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I. INTRODUCTION 



High accuracy measurements of atomic lifetimes are difficult, with only a few examples 
of precisions well under one percent known. For this reason, consideration of radiative cor- 
rections to lifetimes is usually not necessary. Until recently the most prominent exception 
to this situation was the decay rates of orthopositronium and parapositronium, where de- 
terminations of accuracy 180 ppm [lj and 215 ppm [2j respectively have been made. In 
these exotic atoms the radiative corrections start in order a with large coefficients, and the 
leading radiative correction are clearly visible, both because in these two-body atoms the 
wave function is known analytically and because quantum electrodynamic (QED) correc- 
tions start in order a with large coefficients, giving contributions of 2.1 percent and 0.5 
percent respectively. 

In recent years, however, a new approach has been developed that exploits the fact 
that the dipole-dipole potential between two alkali atoms, which goes as as C3/R 3 , can be 
accurately measured, and C3 is proportional to the lifetimes of np states of the atoms. This 
las allowed the determination of the lifetime of the 2pi/ 2 state of lithium as 27.102(7) ns 
jij], the lifetime of the 3p 3 / 2 state of sodium as 16.230(16) ns 0, the lifetimes of the ^pi/ 2 
and 4p 3 / 2 states of potassium as 26.69(5)ns and 26.34(5) ns [5J, the lifetimes of the 5pi/2 
and 5^3/2 states of rubidium as 27.75(8) ns and 26.25(8) ns p, and the lifetimes of the 
Qpi/2 and 6p 3 / 2 states of cesium as 34.88(2) ns and 30.462(3) ns Q. (References to other 
experimental determinations of lifetimes, some of which are of higher accuracy, can be found 
in the above references. We note also that the lifetimes of the francium 7pi/ 2 and 7p 3 / 2 states 
have been measured as 29.45(11) ns and 21.02(11) ns using a different technique |8].) These 
high accuracies, which for matrix elements correspond at the best to 50 ppm, now make the 
calculation of radiative corrections of interest, even though unlike the case of positronium 



the- corrections to El matrix elements arc known to enter, for hydrogenic ions jjj, hi order 
a(Za) 2 . 

There are two other reasons for carrying out calculations of these radiative corrections. 
Firstly, this is a relatively unexplored region of QED. It has only been very recently that 
the first full calculations of one-loop radiative corrections to the decay rate of the hydrogen 
isoelectronic sequence have been carried out jjj. That calculation was done by considering 
the imaginary part of the two- loop Lamb shift, which is equivalent to calculating the shift in 
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the lifetime. Here we adopt a different method, and instead calculate radiative corrections 
to the associated transition matrix element. The introduction of this technique requires a 
nontrivial modification of the previous formalism, and also provides a check of that method. 

The second additional reason for calculating radiative corrections to transition matrix 
elements in the alkalis is the interest in parity n_ation (PNC) in eesitnn Q. Re- 
cently a large binding correction to the Z boson-electron vertex radiative correction has been 
found that has significant implications for the standard model. The lowest-order radiative 

nnn 

correction, —a/2ir, has been shown I12L Il3| to be enhanced through binding corrections 
by almost an order of magnitude. However, the actual radiative correction is to the El 
matrix element of a 6s electron to a 7s electron, with an opposite parity component present 
in one of these electrons induced by Z exchange with the nucleus. The full radiative cor- 
rection calculation needed then involves the evaluation of the diagrams in Fig. ^ of which 
diagram^;, where the radiative correction is on the photon rather than the Z vertex, is of 
the kind to be treated here. Thus the calculations on cesium that will be presented here will 
be of use for this larger scale task. 

We note that the methods used in this paper have previously been applied to calculating 
the Lamb shift 14] and the radiative correction to hyperfine splitting |l5( for the ground 
states of the alkalis, where QED effects larger than experimental accuracy were found. 
However, the calculations remain untested because of the relative inaccuracy of many-body 
methods (with the exception of lithium, where highly accurate variational methods are 
available 16|). We will see that the same situation is present for alkali lifetimes, and thus 
the present work provides yet another impetus to many-body theory to reach the accuracies 
presently of interest for both radiative corrections and experiment. 

In the next section, the generalization of the method used for the previous calculations of 
radiative corrections in the alkalis will be laid out. Of particular interest are certain issues 
related to the fact that the Gell-Mann-Low formalism for energies used in the past work 
has to be changed because we are now instead dealing with matrix elements. In Section ITTTj 
the technique is applied to hydrogenic ions to compare with previous work, and in Section 
TVl the main results of the paper are presented. We conclude with a discussion of how this 
approach can be made more accurate and how it can be generalized for application to the 
calculation of radiative corrections to cesium PNC. 
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II. FORMALISM 



When one is interested in calculating energy shifts in atoms, Sucher's generalization [17 1 
of the Gell-Mann-Low formalism provides a systematic way to derive them from the S-matrix 
through the formula 

A£= I*" -^-^ (1) 
where is the S-matrix with the interaction Hamiltonian Hj(t) multiplied by the factor 

Ae-^L (2) 

The overall factor of e is compensated by the fact that the S-matrix diverges as 1/e. In higher 
orders, factors of 1/e 2 are encountered in the numerator that are canceled by 1/e terms 
coming from expanding the denominator. When considering matrix elements the factor A, 
which accounts for a combinatorial factor specific to the Gell-Mann-Low formalism, will not 
be used, but we continue to use the exponential damping factor e, which will lead to the 
frequent occurrence of what is effectively the delta function, 

D &) = --Am- ( 3 ) 

71 x 2 + e z 

We replace the time-independent perturbation of Ref. [3] with a time dependent Hamil- 
tonian in the length-gauge form appropriate for describing the absorption of a photon in a 
electromagnetic field of strength Eq linearly polarized in the z direction, 

H LG = eE Jd 3 x^ ] (x,t)x-ze- iuJt ip(x,t)a(k 1 z). (4) 

We have made the dipole approximation, a good approximation for the neutral alkalis, so the 
momentum of the initial photon, k, plays no role in the following. In addition we suppress 
the factor eEo in the following. We consider the matrix element of this Hamiltonian between 
an initial state v taken to be a ns ate (2s for lithium, 3s for sodium, etc.) along with a 
photon with energy u, and a final state w taken to be a npi/ 2 or np 3 / 2 state with the same 
principal quantum number n. It is important to keep the photon energy uo distinct from the 
resonance energy ujq = e w — e v , and in particular the limit e — > is always understood to be 
taken before the limit u — > Uq. In lowest order, the S-matrix is then 

S = —2niD t (uj — u )z wv , (5) 
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and following the convention just mentioned, this becomes the usual 



S = [— 2itid(u — uq)]z wv . (6) 

In the following, we shall suppress the factor in square brackets. In Table HI we present 
results for the dipole matrix element r wv , with 



dr r 

o 



9w(r)g v (r) + f w {r)f v {r) , (7) 



from which z wv can be obtained by either multiplying by a factor of (l/3)(— \y™~ mv or 
-\/2/3 for the ns\/2 —npi/2 and nsi/2 — np^/2 transitions, respectively. These are lowest-order 
results obtained using the local Kohn-Sham potential, a description of which can be found 
in Ref. Q- Issues involved in correcting these results with many-body methods will be 
addressed in the concluding section, but here we will concentrate on radiative corrections. 
In evaluating these we will work in terms of the ratio R wv defined through 

(X 

$ z wv — z wv R wv ^ (8) 

7T 

such that 

^wv &Zwv Zwv ( 1 Rwv ) • (9) 

V 7T / 

At this point we also define the frequently occurring self-energy operator 

//* f d n k 1 — 
d 3 xjd 3 yj -j—^ k 2 + iS MxhtiS F {x, y; e - k )Y^j{y), (10) 

in terms of which the lowest-order self-energy part of the Lamb shift of a state v, treated in 
Ref. jl^, is simply S w (e„). 

Most of the discussion of radiative corrections to z wv is very similar to the treatment 
using the Gell-Mann-Low formalism given in Ref. Q], except here we pull out a factor of 
—27ii5(uj — uj ) as opposed to isolating 1/e terms. As described in more detail in that paper, 
three diagrams shown in Fig. El contribute to the radiative correction to the matrix element. 
The S-matrix associated with the vertex (V) diagram of Fig. |2Jd is given by 

x S F (y, z- E 2 ) 1 ^ v {z)D t {E l + k - e w )D e {E 2 + k - e v )D t {E 2 -E l +uo). (11) 
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One can make the substitutions E\ — > e w — k Q and E 2 — > e v — ko in the electron propagators, 
which then allows the E% and E 2 integrations to be carried out, 

fd^fdE2 D ^ + tw)D ^ E2 + ko _ ev ) De ( E2 -E x -co) 

J 111 J ATI 

3e 1 



4tt 3 (u - luq) 2 + 9e 2 
6(u- u ) 



(12) 



471- 2 

We note that in energy calculations where the factor u — ujq vanishes, a factor 1/3 results 
that is canceled because a derivative with respect to the factor A present in the energy 
formula acts on a factor A 3 . Here the factor 1/3 is not present for a different reason, that 
being the fact that the effective infinitesimal factor used to obtain the energy conserving 
delta function is 3e. We then find, after pulling out the factors mentioned above, a vertex 
contribution of 

d n k e&V-*) 
(2ir) n k 2 + i5 

x S F (y, z\ e v - koh^ipviz). (13) 



d 3 x J d 3 y J d 3 z J — — ^ - — ip w (x)~i^S F (x, y; e w - k )y-z^ 



This expression has both ultraviolet divergences and reference state singularities that cancel 
with the side graphs of Figs.Ek and|2t discussed below. The ultraviolet divergence is isolated 
analytically by replacing both bound state propagators with free propagators, which gives an 
ultraviolet divergent term along with a finite remainder we tabulate as R WV (V; 00) in the first 
row of Table ITT1 (We note at this point that we present results only for ns — npi/ 2 transitions. 
Those for ns — np 3 / 2 transitions should not be too different for neutral systems considered 
here, as Native eorreetions to the Methnes of the 2p 1/2 and 2 P3/2 states are essentially 
the same for low-Z hydrogenic ions jSj.) We then form the ultraviolet finite difference of 
Sz wv (V) and 6z wv (V;00) and evaluate it in coordinate space. A Wick rotation k — > iu 
is carried out, which passes poles: the separate contributions are tabulated as R wv (V;iu) 
and R WV (V; Poles) in the second and third rows of Table ITT1 The reference state singularity 
mentioned above is present in the first term, and is regulated through the replacement 
£ v — iS) and e w — > e w (l — id) with 5 typically chosen to be 10~ 6 . 

The treatment of the side graphs also differs from the previous energy approach. The 
starting expression for the "side-right" (SR) diagram of Fig. is 

/f { [ d n k e ^ c dE fdE - 
d 3 xjd 3 yjd 3 zj ^ - w - [ ^ j^j^ US)lo x ■ z S F (x, y; E x )^ 



x S F (y, z; E 2 )YM^) D e{Ei + to - t w )D e {E 2 + k - t v )D e {E 2 + k - E{). (14) 
Replacing the first electron propagator with a spectral representation gives 

C a- 3 V- [j3 [ dE l f dE 2 ^ w (x)-f X-Z1p m (x) 
Ssr = Si^a 2^ / d J x / -— / — — Y, mv (E 2 ) 



m 



2tt J 2tt Ei — e m (l — i5) 
x D e (E x +oj- e w )D e (E 2 + k - e v )D e (E 2 + k - E x ). (15) 

The D e functions emphasize E\ = e v . If the state v is excluded in the sum over states m, the 
same kind of manipulations applied to the vertex graph allow one to determine a "perturbed 
orbital" (PO) contribution from the SR diagram of 

Sz wv (PO;ns) = — Z mv (e v ). (16) 

J ~f e v — e m + to 

The notation (PO; ns) refers to the fact that this contribution is essentially the self-energy 
S^(e„), where v is a perturbation of the ns state. A similar contribution, 5z wv (PO;np) 
arises from the side-left diagram, and is so designated because it is S mii (e m ). Its specific 
value is 

x mn \ [j3 v-v ( ^m{x)~f x -ztfj v {x) 

5z wv (PO;np) = d 6 x ^ S wm (e w ) — — — . (17) 

J m ^ w £ w e m + id 

The case when m = v requires more care. We need to make a Taylor expansion of the 

electron propagator in the self-energy function around the point E 2 = e v — ko, 

S F (y, z; E 2 ) = S F (y, z; e v - k ) + (E 2 - e v + k )S' F (y, z; e v - k ) + ■ • • . (18) 

The first term of the expansion is highly divergent, involving the integral 
f dE p dE 1 

hr 1 hr 1 v D ^ + UJ ~ e ») A(£ 2 + h- e v )D e (E 2 + k - E x ) 

J lit J Zn Ei — e v + to 



1 f UJ — UJq 



(19) 



8m 3 [(u - u ) 2 + e 2 ] 1 7T 3 [(u - u o y + 9e 2 ] [(u - u Q f + e 2 ] j " 
Because we take the limit e — > before to — > loq the second term can be dropped, leaving 
the divergent expression 

i 

bz wv (P>vj) = -— z wv H w {e v ) (20) 

In the energy formalism this corresponds to a 1/e 2 term that is canceled by a term in the 
denominator of Eq. |T|). Here it does not cancel, but instead forms the second term of the 
Taylor expansion of the phase factor 

g— i[E vv (e v )+Y: ww (e w )]/2e (21) 



multiplying z wv , where we have now included the effect of the side- left diagram. Considera- 
tion of higher-order diagrams with two self-energies and a photon interaction show that even 
more divergent terms going as, for example, £^(e„)/e 2 are present that continue the Taylor 
expansion. Thus this divergent term is not directly canceled as in the energy formalism but 
instead can be ignored because it enters only as a phase 19]. However, the second term in 
Eq. (fTHj) does contribute a finite amount to the scattering amplitude, 

6z wv (V) = \ z wv [V m {e v ) + , (22) 

where we have included a similar term arising from the side-left diagram. The treatment 
of these derivative (D) terms follows that of the vertex, but in this case we simply sum the 
finite part of the free propagator term, the Wick rotated term, and the terms in which poles 
are encircled and present the results as R WV (D) in the fourth row of Table ITT1 The ultraviolet 
infinite part of the derivative terms cancels a similar term from the vertex exactly, and the 
reference states are allowed to cancel numerically. Finally, the two perturbed orbital terms 
are listed as R wv (PO;ns) and R wv (PO;np) in the fifth and sixth rows of Table HT1 The 
sums of all contributions give the radiative correction results R wv shown in the last row of 
Table HH 



III. HYDROGENIC IONS 

In a previous paper the imaginary part of the two-loop Lamb shift was used to 
calculate radiative corrections to the 2p 1 / 2 and 2p 3 / 2 lifetimes of hydrogenic ions. In this 
section, we redo the 2pi/ 2 correction for one of these ions using the present matrix element 
formalism as a check of both approaches and to show the role of radiative energy shifts in 
lifetime corrections. The 2p 3 / 2 case included Ml decays, and so cannot be used for a direct 
comparison. 

In Ref. J9J, the radiative correction to the decay rate T is defined in terms of the function 
R{Za) by 

(23) 

An important feature of these corrections is that both radiative corrections to energies and 
radiative corrections to matrix elements are of similar importance. As El decay rates in 
the length gauge are the product of the cube of the energy difference u and the square 
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l + -R(Za) 

71 



of the transition matrix element z wv , one component of R(Za) which arises from energy 
corrections should be given by 



3— = -(Za) 

hi 7T 



32 

-—\n(Za)~ 2 + 22.815 



(24) 



where we have used the standard values of the leading Za-expansion contributions to the 
Is and 2px/2 Lamb shifts. At Z — 5, this contributes —0.064 to R(Za). However, the 
actual value of this function calculated in Ref. jjj] was —0.014, so that a substantial positive 
contribution of about 0.050 must come from the shift in the matrix element. We have 
carried out this calculation using the techniques described above and find a result of 0.025 
for the function R wv defined in Eq. (jHJ). This result is in perfect agreement with the expected 
value, as it has to be doubled to account for the fact that the matrix element enters as a 
square in decay rate calculations. But while this result shows the consistency between the 
present matrix-element and the former decay-rate approaches, there is no improvement in 
the numerical accuracy using our present method for these radiative correction calculations. 



IV. RESULTS AND DISCUSSION 

As found with our previous work on El matrix elements for hydrogenic ions j^, there 
is a high degree of cancellation between contributing terms shown in the first six rows of 
Table [H] This is due to the fact that radiative corrections to decay rates enter in order 
a(Za) 2 . In fact, if one puts the perturbed orbital contributions aside, the cancellation is 
almost complete. As a result, the radiative correction is dominated by the (PO, ns) terms, 
which are the larger of the two perturbed orbital terms. Since it is a standard practice with 
the alkalis to derive matrix elements from observed lifetimes by dividing out the cube of the 
experimental energy differences, we do not need to consider radiative corrections to energies 
here, as they are automatically excluded in empirically extracted matrix elements. 

The experimental accuracies for the alkali matrix elements are half the accuracies of the 
lifetimes quoted in the introduction, specifically being 130 ppm, 490ppm, 940ppm, 1440ppm, 
50 ppm, and 1870 ppm for lithium through francium. Comparing with the results of Table ITTf 
we see corresponding theoretical contributions of 12 ppm, 39ppm, 149ppm, 488ppm, 606ppm, 
and 2158 ppm. Thus these contributions are in principle measurable for several of the alkalis, 
particularly cesium and francium. 
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For the case of lithium, the present approach is certainly less accurate than that avail- 
able from nonrelativistic quantum electrodynamics (NRQED). The nonrelativistic lithium 
wavefunction is extremely well understood using variational techniques and the same 
kind of NRQED calculations reported in Ref. jjj for hydrogenic ions can certainly be used 
for neutral lithium. A nonrelativistic evaluation of the lifetime of the lithium 2»i/ 2 state 

n 

has been carried out in Ref. [20j with a result of 27.1045(14) ns which includes finite mass 
corrections along with relativistic corrections. The quoted error came mainly from the fact 
that the relativistic corrections were not directly calculated, and could be eliminated with 
rigorous relativistic calculations. Regardless, the radiative correction calculated here is so 
small that it should not make any difference to the theoretical lifetime in comparing with 
experiment. 

The cases of cesium and francium are more problematic because of the complexity of their 
wave functions. A great deal of effort has gone into treating these wave function accurately, 
largely spurred by interest in PNC transitions. However, while experimental accuracy is now 
certainly high enough for detection of radiative corrections of the size found here for cesium 
and almost at the level needed for francium, further advances in many-body theory will be 
required, as was the case for the Lamb shift and the hyperfine splitting [ljj, before one 
can decisively say this effect has been seen. Incidentally, a theoretical advantage for lifetimes 
as compared to hyperfine splittings is the relative insensitivity of lifetimes to nuclear effects: 
as emphasized in Ref. uncertainties in the distribution of nuclear magnetism lead to 
theoretical uncertainties that are difficult to control. 

A major spur to theoretical work on accurate cesium wave functions is the fact that the 
observation of PNC in the atom |l(| has significant implications for particle physics. Before 
the large binding correction to the lowest-order radiative correction —a/2Tr was found, a 
discrepancy with experiment existed. As discussed in the introduction, a motivation for the 
present work was to provide a basis for the full calculation of the diagrams shown in Fig. 
We note, for example, that diagram ^ is of the same form as the vertex diagram V in 
Fig. Eh, with the difference that the 6s (or 6pi/ 2 ) state in diagramEb must be replaced with 
a perturbed orbital of the same parity which arises from the 6pi/2 (or 6s) state perturbed 
by a Z boson exchange with the nucleus. Thus a relatively straightforward modification of 
the codes developed for the present calculations will allow the determination of this term. 
In fact, with the exception of diagram le, the entire PNC calculations of Fig. ^ can be 



10 



treated as perturbations of either the weak interaction calculations ^J, the self-energy 
I I 

calculations [14] , or the present transition matrix calculations. We are thus in a position to 
carry out the bulk of the full radiative correction calculations. The fact that the previously 
neglected radiative correction to the El matrix element is 0.261a/7r as shown in Table ITD 
makes it likely that the full PNC calculation will differ quantitatively, though perhaps not 
qualitatively, from the calculations that consider only the Z vertex. 

Since the lowest-order results presented in Table |U are substantially corrected by higher- 
order effects in many-body perturbation theory (MBPT) 2lj], similar corrections can be 
expected for the radiative corrections. The simplest way to account for these corrections 
is to assume that R wv , defined as a ratio of the radiative correction to the lowest-order 
matrix element, is valid when the lowest-order dipole matrix element is replaced with more 
accurate results obtained with MBPT methods. As the size of the radiative corrections 
is small, even if this is only roughly true, the basic size of the effect would have been 
established here. To attain more accuracy, a QED perturbation theory approach could be 
taken. In that case, one would first consider graphs of the type shown in Fig. ^ where the Z 
vertex is replaced by interactions with other electrons. For even higher accuracy, one could 
consider yet more complicated graphs with one absorbed photon, one radiative photon, and 
two interactions with the other electrons: this is known, in the case without the radiative 
photon, to give results within a few percent of the experimental answers 2l|. Progress in 



these large scale calculations, taken together with continuing advances in experiment and 
many-body methods, should allow tests of the lifetimes of the alkalis at the level achieved 
for positronium. 
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TABLE I: Dipole matrix elements r wv for ns 1 / 2 — np l / 2 and ns 1 / 2 — np 3 / 2 transitions in the alkalis 
with Kohn-Sham potentials. Units a.u.. 



transition 


Li 


Na 


K 


Rb 


Cs 


Fr 


ns l/2 - np l/2 


4.171 


4.588 


5.681 


6.009 


6.585 


6.511 


ns 1/2 - np 3/2 


4.171 


4.587 


5.679 


5.995 


6.545 


6.328 


TABLE II: Self-energy contributions R wv 


to El matrix elements for nsi/ 2 


- nPi/2 


transitions in 


the alkalis: error 


of 0.001 for Sum. Units (a/ir)z wv . 








Term 


Li 


Na 


K 


Rb 


Cs 


Fr 


#™(V;00) 


-9.162 


-9.308 


-9.509 


-9.573 


-9.655 


-9.647 


R wv (V;iw) 


134.712 


105.062 


127.516 


125.136 


1 QO 71 Q 


1 1 9 zLfizL 


R WV (V; Poles) 


-137.256 


-107.481 


-129.741 


-127.273 


■134.877 


-114.624 


R WV (D) 


11.698 


11.712 


11.731 


11.738 


11.748 


11.747 


R wv (PO;ns) 


0.003 


0.031 


0.067 


0.182 


0.326 


0.787 


R wv (PO;np) 


0.000 


0.001 


0.000 


0.000 


0.000 


0.202 


Sum 


-0.005 
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(c) (d) 




(e) (f) 



FIG. 1: Feynman diagrams for the radiative correction to electron excitation by a laser photon, 
indicated by the wavy line terminated with a triangle, in the presence of interaction with the 
nucleus through the exchange of a Z boson, indicated by the dashed line terminated with a cross. 
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FIG. 2: Feynman diagrams for the radiative correction to the matrix element for ns + photon 
np. 
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